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ABSTRACT 

Aims. We attempt to determine if a dependency on spectral subtype or vsim exists for stars undergoing phase-changes between B and 
Be states, as well as for those stars exhibiting variability in Ha emission. 

Methods. We analyse the changes in Ha line strength for a sample of 55 Be stars of varying spectral types and luminosity classes 
using five epochs of observations taken over a ten year period between 1998 and 2010. 

Results. We find i) that the typical timescale between which full phase transitions occur is most likely of the order of centuries, 
although no dependency on spectral subtype or vsim could be determined due to the low frequency of phase-changing events observed 
in our sample, ii) that stars with earlier spectral types and larger values of vsini show a greater degree of variability in Ha emission 
over the timescales probed in this study, and iii) a trend of increasing variability between the shortest and longest baselines for stars 
of later spectral types and with smaller values of vsini. 
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The e mission lines d isplayed by Be stars are known to be tran- 
sient (Kogtirel [l990T) . with emiss ion varying over timescales 
from less than a day dPercv et al.l I2002I) to 
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dMiroshnichenko et al.l |2002|) . This emission v ariability is re- 
flected in the first working definition given bv ICoUinsI (I1987I) 
who defined a Be star as "a non-supergiant B star whose spec- 
trum has, or had at some time, one or more Balmer lines in emis- 
sion". 

More recently, this definition has had to be refined to ex- 
clude emission-line stars whose disk formation mechanisms 
were accretion, e.g. Herbig AeBe stars. Consequently, a new 
subgroup of these emission-line stars, termed "classical Be 
stars", was created. Classical Be stars are characterised by thei r 
rapid rotation and dec retion disks (see iPort er & Rivinius 2003I) . 
Martavan et al.l d201 ll) proposed a refinement of Collins' defini- 
tion to "... a star with innate or acquired very fast rotation which 
combined to other mechanisms . . . leads to episodic matter ejec- 
tions creating a circumstellar decretion disk or envelope" to not 
only reflect the larger range of spectral t ypes over which the 
phenomena have been observed to occur (Conti & Leedll974 



Rauw et al. 2007) albeit with lesser frequency (Negueruel a et al.l 



2004), but also to constrain the method whereby the disk is 



formed. It is these classical Be stars that are the subject of the 
following study. 

In 1999, a multiwavelengt h survey of a representative sample 
of 58 Be stars was undertaken ([Steele et all 19991: IClark & Steele] 
2000: ISteele & Clarkll200lHHowells et allboOll) . The sample is 
"representative" in that it attempted to include objects of each 
spectral type and luminosity class and so doesn't serve to reflect 
the spectral type or luminosity class distribution of Be stars. 

One of the conclusions drawn from this survey was that at the 
time of observation, non-emission line stars in the sample (i.e. by 
definition, ones that must have shown emission in the past) had 



lower rotational velocities and earlier spectral types, an implica- 
tion being that they m ay be more prone to B/Be phase changes 
jSteele & Clarkl 120011) . Using optical spectrograph^ follow up 
taken ten years after the sample was first observed, the validity 
of these statements will be investigated further in this paper. 



2. Observations and data reduction 



rmb@astro.livjm.ac.uk 



The targets observed were taken from Stee le et aD (Il999l here- 
after S99). The catalogue names of the targets and some of their 
fundamental properties (viz. spectral class, luminosity class, V 
band mag., vsim) are shown in TableQ] 

In S99, the stars were r eclassified under the MK 
scheme ( Morg an" & Keenarlll973l) using spectra taken by the 
Intermediate Dispersion Spectrograph (IDS) on the Isaac 
Newton Telescope (INT), rather than relying on values taken 
from the literature. Rotational velocities were also determined in 
S99 by fitting Gaussian profiles to 4 Hel lines at 4026A, 4143A, 
4387A and 4471 A and applying the full w idth half maximum 
- vsin; correlation of Slette bak et al.l (119751) . It should be noted 
that this method of inferring a vsini has not taken into account the 
effect of gravity darkening, which has been shown to introduce 
a redundancy b etween line profile width and vsim at the largest 
rotation speeds dTownsend et al. 2004; Fre mat et al.l20 05), lead- 
ing to underestimates of the true rotation speed for the fastest 
rotators. Residual emission within the Hel lines may have also 
introduced a bias towards lower vsini for earlier subtypes. 

Five epochs of observations were used in the follow- 
ing analysis (hereafter referred to as the INT, Fl, F2, F3 
and F4 datase ts) and were obtained using both the IDS and 
FRODOSp ec dMorales-Rueda et alj [2004) on the Liverpool 
Telescope (ISteele et al.ll2004l) . 

The IDS observations used in this analysis were made on 
the night of 1998 August 3 using the R1200Y grating with a 
slit width of 1.15 arcsec, corresponding to a dispersion of ~0.5 
A/pixel on the EEV12 CCD. A central wavelength of 6560 
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Table 1. The sample targets with their corresponding stellar parameters taken from S99. "sh" indicates that shell lines are observed 
in the spectra. The measured equivalent widths of the sample Ha lines are shown for each epoch. Missing observations are marked 
with an ellipsis. 



Object Name 


Spec. Type Lum. Class 


VBand 


vsim 
(km/s) 


EW 
INT 


Fl 


F2 


F3 


F4 


CD -28 14778 


R? 

DZ. 


TTT 

hi 


8.95 


153 


-32.9 






-31.0 




CD -27 11872 


Rfl S 

DUJ 


V TTT 

V-lll 


8.69 


224 




-34.9 




-37.5 




CD -27 16010 


R8 

Do 


TV 

1 V 


4.20 


187 


0.0 


-1.6 


-1.0 


-2.0 


-1.1 


CD -25 12642 


Rfl 7 


TTT 
in 


9.00 


77 


2.7 


2.5 




2.6 




CD -22 13183 


R7 


V 

V 


7.90 


174 


-8.3 


-14.2 




-14.5 




BD -20 05381 


RS 

D J 


V 

V 


7.80 


202 


-12.1 


-13.3 




-14.6 




BD -19 05036 


RzL 
D + 


TTT 
111 


7.92 


121 


4.5 


4.6 








BD -02 05328 


R7 

D / 


v 

V 


6.22 


151 


-6.3 


-5.2 




-5.8 


-6.3 


BD -01 03834 


R? 


TV 

1 V 


8.14 


168 


-35.9 


-50.1 




-45.2 




BD -00 03543 


R7 

D / 


V 
v 


6.88 


271 


-1.8 


-5.0 




-6.4 




BD +02 03815 


R7 S 


SI1 


6.92 


224 


-5.5 


-6.2 




-7.4 




BD +05 03704 


R? S 


V 

V 


6.13 


221 


0.3 


4.0 




3.8 




BD +17 04087 


Rfi 

DU 


TTT V 

111- V 


10.6 


156 


3.9 


4.2 




3.4 




BD +19 00578 


RR 

Do 


V 

V 


5.69 


240 


4.0 


5.4 


5.7 


4.9 


4.6 


BD +20 04449 


Rf) 

DU 


TTT 

in 


8.3 


81 


2.4 


4.0 


3.4 


3.2 


3.0 


BD +21 04695 


Rfi 

DU 


TTT V 

111- V 


5.78 


146 


-0.7 




-14.2 


-15.3 




BD +23 01148 


R? 


TTT 

in 


7.36 


101 


0.6 




0.7 




1.1 


BD +25 04083 


Rf) 7 

DU. / 


TTT 
in 


8.94 


79 


2.5 


2.5 


3.3 


2.5 


2.4 


BD +27 00797 


Rfl S 

DU. J 


V 

V 


9.86 


148 


-18.2 






-29.7 


-22.7 


BD +27 00850 


R1 S 

D 1 . J 


TV 

1 V 


9.38 


112 


4.0 




4.2 


3.1 


2.7 


BD +29 03842 


R1 

D 1 


TT 
11 


10.11 


91 


1.6 


2.4 




1.7 


1.7 


BD +29 04453 


R1 S 

D 1 . J 


v 

V 


8.10 


317 


-54.3 


-61.7 


-60.1 


-58.8 


-63.8 


BD +30 03227 


RA 

D4 


V 

V 


6.58 


218 


1.8 


5.6 


5.6 


6.0 




BD +31 04018 


R1 S 

D L.J 


V 

V 


7.16 


211 


-18.6 


-21.1 




-19.6 


-20.9 


BD +36 03946 


R 1 

D 1 


V 

V 


9.20 


186 


-23.2 


-31.7 




-35.7 


-35.6 


BD +37 00675 


R7 

D / 


V 

V 


6.16 


207 


-12.7 


-14.0 


-12.2 


-11.0 


-11.2 


BD +37 03856 


Rfl S 

DU. J 


V 

V 


10.20 


104 


3.6 


3.4 




3.3 


2.7 


BD +40 01213 


R? S 


TV 

1 V 


7.34 


128 


-12.2 




-14.6 


-17.9 


-17.3 


BD +42 01376 


B2 


v 


7.28 


196 


-7.8 


-15.2 


-15.9 


-12.8 


-13.9 


BD +42 04538 


g2.5 


v 


8.02 


282 


-37.5 


-37.8 


-38.8 


-41.6 


-41.4 


BD +43 01048 


Rft 

DU 


TTT 

lll-Sll 


9.53 


220 


-6.0 




-8.8 


-10.3 


-9.5 


BD +45 00933 


B1.5 


V 

V 


8.06 


148 


3.9 


3.5 


3.7 


3.7 




BD +45 03879 


R1 S 

D 1 . J 


V 

V 


8.48 


193 


-31.1 


-23.4 


-20.6 


-22.7 


-25.0 


BD +46 00275 


R*i 

D J 


III 


4.25 


113 


1.0 


0.9 


0.8 


1.5 


2.4 


BD +47 00183 


R? S 


V 

V 


4.50 


173 


-29.7 


-36.1 




-38.0 


-37.2 


BD +47 00857 


34 


V-IV 


4.23 


212 


-36.3 


-38.7 


-38.6 


-43.0 


-41.1 


BD +47 00939 


R? s 


v 

V 


4.04 


163 


-22.4 


-25.6 


-24.5 


-25.4 


-26.1 


BD +47 03985 


B1.5 


sh 


5.42 


284 


-17.8 




-23.2 


-25.3 




BD +49 00614 


RS 

D J 


TTT 

in 


7.57 


90 


-4.9 




-6.4 


-7.5 


-7.0 


BD +50 00825 


R7 

D / 


V 

V 


6.15 


187 


-1.5 


-2.1 


-1.7 


-2.6 


-2.1 


BD +50 03430 


RR 
Do 


v 

V 


7.02 


230 


-6.2 


-9.8 


-9.7 


-10.6 


-10.9 


BD +51 03091 


B7 


III 


6.19 


106 


2.2 


1.0 




0.7 


1.1 


BD +53 02599 


B8 


V 


8.08 


191 


-0.7 


1.2 


1.8 


1.8 


1.8 


BD +55 00552 


B4 


V 


7.90 


292 


-1.3 


-9.4 


-10.4 


-10.8 


-10.6 


BD +55 00605 


Bl 


V 


9.34 


126 


-5.9 


-7.5 


-5.2 


2.2 


-0.7 


BD +55 02411 


B8.5 


V 


5.89 


159 


2.9 


2.1 


2.7 


2.9 


2.7 


BD +56 00473 


Bl 


v-m 


9.08 


238 


-37.4 


-18.2 


-20.1 




-17.3 


BD +56 00478 


B1.5 


V 


8.51 


157 


-9.0 


-12.5 


-13.5 


-16.8 




BD +56 00484 


Bl 


V 


9.62 


173 


-47.8 


-46.6 


-47.2 


-44.6 


-36.7 


BD +56 00493 


Bl 


V-IV 


9.62 


270 


2.0 


-0.3 


1.0 


-2.4 


-1.9 


BD +56 00511 


Bl 


III 


9.11 


99 


-5.1 


-4.2 


-5.4 


-6.4 


-6.1 


BD +56 00573 


B1.5 


V 


9.66 


250 


-67.5 


-38.5 


-73.3 


-77.1 


-74.6 


BD +57 00681 


B0.5 


V 


8.66 


147 


0.2 


1.4 


1.6 


0.6 




BD +58 00554 


B7 


V 


9.16 


229 


-9.5 


-7.9 


-7.6 


-8.9 




BD +58 02320 


B2 


V 


9.51 


243 


-4.1 


-3.9 


-5.2 


-5.9 


-6.0 
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A was chosen, giving a wavelength range between 5800 - 7100 
A. 

Reduction of INT data was performed using Figaro accord- 
ing to the standard prescription for long-slit spectra. Tungsten 
lamp flats were used to map the slit response and sky flats the 
pixel-to-pixel flat field variations. Target spectra were extracted 
using simple extraction and a sky region of the same size was 
also extracted and subtracted from each target spectrum. To 
wavelength calibrate the target frames, arc exposures were taken 
at the start, middle and end of the night using CuAr and CuNe 
lamps. After extracting their spectra, the calibrations were deter- 
mined and copied onto the target spectra. 

All FRODOSpec observations were taken over a period be- 
tween September 2009 and November 2010. Observations for 
the Fl dataset were made between 08/09/2009 and 26/09/2009, 
F2 between 03/1 1/2009 and 26/12/2009,F3 between 18/07/2010 
and 29/08/2010 and F4 between 10/10/2010 and 23/1 1/2010. 

The first two epochs of data, Fl and F2, were taken using 
the red diffraction gratings. The last two, F3 and F4, were taken 
using the higher resolution red VPH gratings. FRODOSpec pro- 
vides wavelength coverage from 3900-5700A (blue arm) and 
5800-9400A (red arm) for the lower resolution configuration 
(R = 2600 and 2200 respectively) and 3900-5 100 A (blue arm) 
and 5900-8000A (red arm) for the higher resolution (R = 5500 
and 5300 respectively). FRO DOSpec reduction was performed 
using the pipeline discussed in Barnsle v et al.l d2012l) . 



certainty, sets of repeated observations of B stars with different 
brightnesses spanning the sample range (V = 5, 8.5 and 10) were 
taken. These observations were taken sequentially on the same 
night with both FRODOSpec's resolution configurations, using 
corresponding exposure times similar to those used to observe 
the sample. 

The standard deviation of the measured Ha absorption EWs 
was calculated for each set of target spectra. The V band mags, 
and the measurement errors for each target are shown in Table [2] 

These errors were associated to all measurements of EW, 
matching them to the nearest target brightness. In the absence 
of equivalent INT observations, IDS spectra were assigned the 
same error figures as those used for FRODOSpec. In analysis 
requiring the photospheric contribution to be subtracted off, the 
errors in measuring the EW and the error in assigning a photo- 
spheric EW were added in quadrature. 

Table 2. V band mags, and values of EW used as the continuum 
normalisation error. 



Object Name 


V 


EW(Hcy) o- 


HR5834 


5.00 


0.30 


HD 75750 


8.48 


0.1 


HD 100340 


10.12 


0.1 



3. Method 

Emission in the Balmer series is commonly used to provide an 
insight into the circumstell ar environment surrounding Be stars 
(Grundstrom& Gies 2006) and results from recombination of 
photospheric radiation within the disk. In this paper, only the 
line strength of Ha is measur ed. Ha is observed i n Be stars with 
a variety of profile shapes dBaneriee et al . 2000), ranging from 
single and double peaked emission to absorption, where no iden- 
tifiable emission component can be detected (i.e. the star is not 
currently in an emission phase). 

3.1. Quantifying the extent of disk loss/formation 

Due to in-filling of photospheric absorption lines, visual inspec- 
tion alone is not reliable at being able to discern between a star 
that shows weak emission and a star showing no emission at 
all, especially when the observations vary in resolution. It was 
therefore necessary to use a more quantitative method, capable 
of quantifying the extent of disk emission in a more robust man- 
ner. For this reason, the equivalent width (EW) of the line was 
used to both assess phase change and quantify the variability of 
the emission. 

In order to measure the EWs, the spectra were first nor- 
malised by fitting a polynomial to the shape of the underlying 
continuum. The EW was then calculated over a suitable wave- 
length range encompassing the line. The EW measurement was 
performed in IDL using the line_eqwidth routine from the 
ltools packag4_t 

The dominant source of error in measuring the EW is the 
determi nation of the con tinuum level and subsequent normal- 
isation dJones et al. I l201lh . with any scaling error in determin- 
ing the level of the continuum leading directly to a multiplica- 
tive error in the EW. To obtain an empirical estimate of this un- 

1 http://fuse.pha.jhu.edu/analysis/fuse_idl_tools.html 



Before EW measurements of the Ha line could be used to 
assess if an object exhibited complete disk loss/formation, it was 
required that the photospheric contribution to the EW from the 
central star, W p kou be subtracted off the measurement of the total 
EW of the line, W, oU i.e. 

W disk = W tot - W phot (1) 

where Wdisk 1S the equivalent width of the emission line arising 
from the disk only. 

W p hot can be established by either using models of stellar 
atmospheres (e.g. lKuruczl[l979l) or by measuring the EW for 
a set of B stars (with no history of emission) over a range of 
spectral types and luminosity classes and assuming a linear fit 
to EW and spectral type for each luminosity class. In order to 
retain model independency in the following analysis, the latter 
method was prefered and a set of B stars, shown in Table [3] were 
observed. 

A plot of the EW against spectral subtype for each luminos- 
ity class is shown in Fig.Q] The corresponding Pearson correla- 
tion coefficients for luminosity classes III, IV and V are -0.97, 
-0.94 and -0.88 respectively, indicating that the data are reason- 
ably well described by a linear fit. 

As an aside, the measured EWs were compared with those 
determined from using Kurucz models. The model was found 
to underestimate for earlier spectral types and overestimate for 
later, possibly due to the effects of line blending at the 20A res- 
olution of the grid around the Ha line. 



4. Results and Discussion 

The measured equivalent widths of the sample Ha lines are 
shown for each epoch in Table Q] 
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Table 3. The observed standards with their corresponding stellar 
parameters and measured Ha equivalent widths. 

Object Name Spec. Type Lum. Class EW 



HD 886 


B2 


IV 


3.8 


HD 20365 


B3 


IV 


4.9 


HD 22951 


B0.5 


V 


4.2 


HD 23180 


Bl 


III 


3.5 


HD 23288 


B7 


IV 


6.1 


HD 23338 


B6 


V 


5.7 


HD 23850 


B8 


III 


5.8 


HD 29248 


B2 


III 


3.4 


HD 147394 


B5 


IV 


5.3 


HD 160762 


B3 


IV 


4.1 


HD 176437 


B9 


III 


6.9 


HD 184930 


B5 


III 


4.9 


HD 186882 


B9.5 


III 


7.6 


HD 196867 


B9 


V 


7.6 


HD 207330 


B2.5 


III 


3.9 


HD 214923 


B8 


V 


7.3 


HD 214993 


B1.5 


III 


4.0 


HD 218376 


B0.5 


III 


3.0 


HD 219688 


B5 


V 


4.5 


HD 222661 


B9.5 


V 


9.6 



1 2 ~~ 
1 1 - 



10 




4 

B Spec 



6 

Subtype 



Fig. 1. Photospheric EW as a function of spectral subtype for lu- 
minosity class III (circle), IV (triangle) and V (square) with lin- 
ear best-fit lines overplotted (solid, dotted and dashed for classes 
III, IV and V respectively). 



-2 



-6 



-10 L 




-10 



-8 



-6 -4-2 2 
Exophotospheric INT EW(Halpha) 



Fig. 2. Determined exophotospheric EWs plotted for the INT- 
Fl baseline. The scale has been selected to show all targets that 
exhibited either disk loss or formation. Targets that have been 
identified as losing their disk are illustrated with filled circles. 
The stability line (y=x) and disk loss/formation (y=0, x=0) axes 
are shown. 



Candidates for disk loss were selected by identification of those 
targets that lay off the stability line (where no EW change could 
reliably be inferred), but whose EW and associated error lay 
within the disk loss (y=0) axis. Likewise, candidates for disk 
formation were selected using similar criteria, except only those 
that lay within the disk formation (x=0) axis were selected. The 
resulting phase-changers are shown with filled symbols in Fig. 
[2] As phase-changing was seen to occur only over the timescale 
of 10 years, only this baseline is shown. 

One of the aims of this analysis was to assess if there was any 
correlation between spectral type and frequency of B/Be phase 
change. Because of the low number of targets that were identi- 
fied undergoing this change, no attempt has been made to de- 
termine if a dependency on this fundamental parameter exists. 
Indeed, complete phase change could only be confirmed at the 
3cr level for two targets in the sample, BD +05 03704 and BD 
+30 03227, with the longest baseline of ten years (INT - Fl). 

If these targets are believed to have truly changed state be- 
tween B and Be phases, it is possible to estimate the approxi- 
mate time over which complete phase change is expected to oc- 
cur. As little observational evidence is available regarding both 
the distribution of phase-changing timescales and the fraction 
of Be stars which undergo complete phase change (most likely 
due to observational bias from incomplete coverage of the time 
domain), the following assumptions have been made: 



4.1. Complete disk loss/formation 

To determine complete disk loss/formation using the calculated 
exophotospheric EWs, values for consecutive pairs of epochs 
(with corresponding baselines of INT - Fl (~10 years), Fl - 
F2 (~2 months), F2 - F3 (~9 months) and F3 - F4 (~3 months)) 
were plotted against each other; see Fig. |2]for an example plot. 



- The time between complete phase changes is randomly dis- 
tributed (neither shorter nor longer times are preferred), and 

- All Be stars are subject to these episodic events of complete 
disk loss/formation. 

With these in mind, a Monte Carlo simulation was set up 
to predict the number of observed phase-changers over a ten 
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year period for different timescales (whilst constraining the min- 
imum time to 10 years - the only baseline over which phase 
change was observed in this sample). The simulation was run 
iteratively using a 55 star dataset, recording how many were ob- 
served to change phase an odd number of times within a defined 
ten year window (an even number would imply that no phase 
change would be observed). In a sample of 55 objects, the max- 
imum time for phase-changing to yield 2 targets observed in a 
ten year baseline was found to be most likely between 500 and 
1000 years. 

However, this maximum time comes with a possible caveat. 
Since Be phenomena are time variable, selection of candi- 
date Be stars for the sample was based on heterogenous data 
taken from a multitide of different sources and observers that 
has been compiled ove r decades to produce a single catalogue 
dJaschek & Egrellll982h . If the original classification of the tar- 
gets used in the sample is untrustworthy, the sample may be 
contaminated by B stars that have never truly shown reliable ex- 
ophotospheric emission. This is certainly a real possibility for 
the fainter targets recorded in the catalogue, and to some extent 
may even be true for some of the brighter targets. 

If fewer true Be targets were present in the sample, the ob- 
served fraction of phase-changers to total targets increases, in- 
creasing the frequency of phase-changes over the timescales sur- 
veyed by this study and effectively decreasing the maximum 
time over which the targets are expected to change phase. To de- 
termine the lower limit for the maximum time, the sample was 
filtered by-eye to remove targets that were always in absorption. 
As noted previously, due to in-filling of photospheric absorption 
lines this is not necessarily indicative that the star is not in emis- 
sion, but serves as an upper limit to the number of stars that are 
non-variable, and thus could have been misclassified. The pre- 
vious Monte Carlo simulation was repeated using the smaller 
dataset, and in this smaller sample of 44 objects, the maximum 
time for phase-changing to yield 2 targets observed in a ten year 
baseline was found to be most likely between 250 and 500 years. 
Finally, we note that this value will drop further if all Be stars are 
not subject to these episodic events of complete disk/loss forma- 
tion. 

4.2. Partial disk variability 

In the following analysis, partial disk variability is conserva- 
tively defined as a 5<x change in EW measurement, where the 
previous constraints placed to force the targets to lie upon either 
the axis of disk loss or formation are removed. As the underlying 
photospheric contribution from the central star will not change 
between epochs, the total EW is used in the analysis of partial 
disk variability. 

4.2.1 . As a function of spectral subtype 

An initial approach to the analysis of partial variability is to as- 
sess its presence over all of the possible baselines. Specifically, if 
an object is seen to vary in at least one of the baselines, it should 
be classed as variable. A plot of the variability over any baseline 
as a function of spectral subtype is shown in Fig. [3] 

Analysis regard ing variability as a func tion of spectral sub- 
type is presented bv lHubert & Floquetl(ll998l) who found that in 
a homogenously obtained Hipparcos sample of 273 Be stars, > 
98% of early Be stars (B0 - B3e) showed some degree of vari- 
ability, compared with only 45% of late Be stars (B7 - B9e) over 
the four years during which the survey was conducted. Using the 
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Fig. 3. Targets showing partial variability over any of the base- 
lines. The hollow histogram represents the distribution of targets 
as a function of spectral subtype. The hatched histogram shows 
the distribution of those targets displaying partial variability. The 
cumulative fractions have been overplotted for those targets dis- 
playing partial variability (solid line) and those not (dashed line). 



same definition of early and late types, similar figures are found 
in this sample, where 84% of the early Be stars in the sample 
were variable with only 46% of the later t ype showing some 
variation at the 5<x level over all baselines. iJones et al.l (|201 ll) 
find the same trend in their sample, but find 45% of their early 
and 29% of late types to be variable. These smaller fractions 
could be the result of the statistic they use to quantify variabil- 
ity, which they state as being "sensitive to significant changes in 
disk density". 

There are several possible explanations for the observation 
of increased variability in earlier spectral types. It could be the 
case that a false-positive misclassification has introduced a bias, 
but such a bias would only alter the trend if later type non- 
variable B stars were preferentially misclassified as being vari- 
able. Filtering the sample by-eye to remove targets that were al- 
ways in absorption (as in 94. U . it is found that the trend remains 
with 100% of the earlier type and 53% of the later type varying. 
Non-variability was found for 7 out of the remaining 45 stars, 
with all 7 stars of spectral type B5 or later. 

Discounting this, the most obvious explanation would be 
that later spectral types are genuinely less variable than earlier. 
However, this may not necessarily be the complete picture, in 
so much that the timescales over which later types vary may not 
have been adequately probed by the timescales in this study. 

To assess if there is any trend in the degree of variability 
across different timescales, pairs of epochs were grouped to- 
gether to construct sets of similar baseline timescales. For the 
purposes of the following analysis, the baselines have been sep- 
arated into < 1 year (F1-F2, F3-F4), about one year (F1-F3, F2- 
F3, F1-F4, F2-F4) and ten years (INT-F1, INT-F2, INT-F3, INT- 
F4). The previous analysis was then redone using the each of the 
sets; the result is shown in Fig. |4] The following figures are pre- 
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sented for both the full dataset and the previous emission-only 
filtered dataset (given in brackets after the full dataset). 
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Fig. 4. Targets showing partial variability over < 1 yr (top), ~ 1 yr 
(middle) and ~10 yrs (bottom). The hollow histograms represent 
the distribution of targets as a function of spectral subtype. The 
hatched histograms show the distribution of those targets dis- 
playing partial variability. The cumulative fractions have been 
overplotted for those targets displaying partial variability (solid 
line) and those not (dashed line). 



For partial variability over timescales less than a year, no 
later types could be confirmed at the 5<x level to be partially 
varying in either dataset, whilst 64% (79%) of the earlier type 
were. For timescales of around 10 years, the number of later 
types varying were found to be 46% in both datasets, with the 
number of earlier types also rising to 84% (100%). This would 
seem to imply some dependency of partial variability as a func- 
tion of spectral type on baseline, but there is little difference in 
the distributions between timescales of one year and timescales 
of 10 years. It is worth noting that as th e distribution of Be stars 
peaks at around a spectral type of B2 (Martavan et al.ll2010l) and 
tails off sharply towards later spectral types, quantitative statis- 
tics derived using these later objects are less certain. 

A two sample Kolmogorov-Smirnov (K-S) test was also used 
to determine if the observed distribution of targets showing par- 
tial variability differed from the underlying parent distribution 
of those which were non-variable as a function of spectral type. 
The calculated rejection probability for the same distribution hy- 
pothesis was found to be 75% (99%), the implication being that 



the stars displaying partial variability may have a different par- 
ent distribution, and thus dependency on spectral type, than those 
which don't. 

4.2.2. As a function of rotational velocity 

A plot of the variability over any baseline as a function of vsim 
is shown in Fig. [5] From visual inspection of the plot, it would 
appear that stars with lower rotational velocity are less variabile 
over the timescales probed by this study. 
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Fig. 5. Targets showing partial variability over any of the base- 
lines. The hollow histogram represents the distribution of targets 
as a function of vsim. The hatched histogram shows the distribu- 
tion of those targets displaying partial variability. The cumula- 
tive fractions have been overplotted for those targets displaying 
partial variability (solid line) and those not (dashed line). 



The two sample K-S test was also used to determine if the 
observed distribution of targets showing partial variability dif- 
fered from the underlying parent distribution of those which 
were non-variable as a function of vsim. The calculated rejec- 
tion probability for the same distribution hypothesis was found 
to be >99% (81%), the implication being that the stars display- 
ing partial variability may have a different parent distribution, 
and thus dependency on vsim, than those which don't. 

For completeness, the trend of variability as a function of 
vsim is also presented using the same epoch pairings discussed in 
34.2.11 the result is shown in Fig. [6] Partial variability is seen to 
occur for nearly 100% of those targets with the highest vsim over 
all baselines, with a trend of increasing variability for longer 
baselines seen for stars with lower vsim. 

However, it is noted that analysis using vsim may be mis- 
leading. To compare stars using this quantity is to compare stars 
which may have different inclination angles (the angle between 
the rotation axis and the line of sight), spectral types and lu- 
minosity classes. A more useful parameter is the ratio of the 
observed rotational velocity to the critical, or break-up, veloc- 
ity of the star (which in itself is a function of spectral type and 
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- The distributions of variable and non-variable stars with re- 
spect to both spectral type and vsin/ are found to be different 
under a two sample K-S test. 
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Fig. 6. Targets showing partial variability over < 1 yr (top), ~ 1 yr 
(middle) and ~10 yrs (bottom). The hollow histograms represent 
the distribution of targets as a function of vsin/. The hatched his- 
tograms show the distribution of those targets displaying partial 
variability. The cumulative fractions have been overplotted for 
those targets displaying partial variability (solid line) and those 
not (dashed line). 



luminosity class) given by to = v rot /v cri ,. However, a model- 
independent determination of v er;( is not attainable, and as such 
we do not attempt to perform analysis using this quantity. 

5. Conclusions 

- The observed frequency at which transitions between the B 
and Be states occur is too low (2/55) for the purposes of as- 
certaining any dependency on spectral type or vsin/. Given 
this, we predict that the typical timescale between which 
these complete phase transitions occur is likely to be of the 
order of centuries. 

- Using a measure of spectral type dependency akin to that 
used by previous authors, we find evidence to suggest that 
stars of an earlier subtype may be more variable. We also 
find a rise in variability for later types (B7-9) when observed 
over baselines greater than a year. However, the weight with 
which this conclusion is made is limited by the small number 
of stars observed at later subtypes. 

- Similar to the previous conclusion, we find that stars with 
larger values of vsin/ may be more variable. Again we note a 
rise in variability for stars with lower values when observed 
over baselines greater than a year. 
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